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Abstract 

The dependence of the one dimensional Bose-Einstein correlation Rid derived from 
proton-proton collisions is examined in terms of energy, hadron multiplicity and number 
of emerging hadronic jets. It is argued that the observed rise of -Rid with energy is the 
result of the increase of the average number of hadronic jets. It is further shown that by 
the relative straightforward measurements of Rib and its chaoticity parameter A one 
should be able to estimate the average number of outgoing hadronic jets emitted over 
the entire final state phase space without the necessity to rely on a particular jet finding 
algorithm and its associated parameters. 
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1 Introduction 



The Bose-Einstein Correlation (BEC) of identical boson pairs, mainly tt^tt^ mesons, produced 
in particles and heavy ion collisions has been studied since the early 1960's [1]. In the 
one-dimensional BEC analyzes one frequently utilizes the Lorentz invariant variable which 
is defined as = —{pi — 'P2Y where pi is the four momentum of the ith identical boson. 
Taking the emitting source of the bosons to be described by a spherical symmetric Gaussian 
density distribution, the correlation function of two identical bosons can then be given e.g. by 

C{Q) = N{1 + Ae-^?DQ^)(l + 5Q + eg2) , (1) 

where Rid (also referred to at times as Rq or i?inv) is the one dimensional (ID) value of the 
particles' emitter source and A is the strength of the BEC effect. This A, which is the 
chaoticity parameter, lies within the range of < A < 1 [2,3]. The parameters 5 and e are 
frequently introduced to account for possible long range correlations, like those related to 
momentum and energy conservations, while is a normalization factor. The results obtained 
from BEC analyzes depend among others on the assumption that the low percentage of the 
non pion outgoing charged tracks does not influence the BEC results. Furthermore the 
particular choice taken from the various reference samples against which the BEC effect is 
compared to has to be kept in mind when common underlying features are searched for in the 
BEC analyzes results. In addition the reported i?iD values for charged identical bosons were 
not always a subject to a Coulomb correction which however is small and estimated to be 
< 0.05 fm ( see e.g. Ref. [4]). In more recent studies the BEC analyzes were extended to two 
and three dimensions to allow for non-spherical emitting sources [2]. 

Our present study concentrates on the ID case which allows us to include early BEC analyzes 
of pp — )■ hadrons collisions carried out at relatively low energies. The main aim of the present 
work is to study the Rid dependence on the pp colliding center of mass energy y^, the 
outgoing charged hadron multiplicity and hadronic jets and to explore the possibility to use it 
for the estimation of the average number of hadronic jets emitted over the full final state 
phase space. 



2 The i^iD dependence on the pp coUiding energy 

The recent operation of the Large Hadron pp Collider (LHC) at CERN has opened up the 
opportunity to examine the Rm derived from BEC of identical pion pairs produced in pp 
collisions at very high center of mass energies up to 7 TeV . A set of Rid values, obtained via 
Eq. (1) and taken from Refs [4-9] are presented in Fig. 1 as a function of the pp colliding 
energy in the range of ^ from 5 GeV to 15 TeV . The errors shown are the statistical and 
systematic ones added in quadrature. From the two ALICE values at 0.9 TeV reported in 
Ref. [4], the one with the considerably smaller systematic error was chosen. In the figure the 
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Figure 1: A set of iiio values and their errors obtained from BEC analyzes, using Eq. (1), of identical 
charged pion pairs produced in pp collisions at several center of mass energies. [4-9] averaged over all 
charge multiplicities. The middle line represents i?iD = 1.64 + 0.14 x ln{^/s) fm as obtained from a 
fit to the data. The dashed lines define the ±1 s.d. error band associated with the fitted parameters. 

CMS and ALICE values at 0.9 TeV are for clarity shifted slightly apart. Here it should be 
noted that the data selection criteria used for the Rid measurements at the LHC energies are 
somewhat different from those used in the BEC analyzes at the lower energies. In particular 
at the LHC energies the observed hadron tracks used for the correlation analyzes are confined 
to a restricted phase space domain commonly specified in terms of the experimentally 
accessible pseudorapidity r) range where r) is defined as 

T] = 0.5 X ln[(|p1 -pl)]. 

Here \p\ is the absolute value of the outgoing particle three momentum vector and is its 
longitudinal momentum component. However it has already been previously observed [10] 
that the BEC lengths show, if at all, only a weak dependence on rj. 

A clear rise of with y/s is seen in Fig. 1 and is here represented by the continuous line 
which was obtained from a fit of the expression 

= a + bx ln{y/s) fm (2) 
2 



to the data yielding a = 1.64 ± 0.11 fm and b = 0.14 ± 0.02 fm with a /dof of 0.3/5 where 
■\/s is given in TeV units. The deviation of ± 1 s.d. is indicated by the dashed hnes in the 
same figure. At first one is lead to presume that the -Rid increase is a direct consequence of 
the energy rise however, this apparently cannot be the prime reason as has been recently 
clearly demonstrated by the CMS collaboration [9]. In their BEC study they have compared 
the i?iD values measured in pp collisions at the two energies of 0.9 and 7.0 TeV as a function 
of the outgoing charge multiplicity A'^chwhich are shown in Fig. 2. The dotted and solid curves 
in the figure are the results from the fits of the relation [11] -RiD(A^ch) = k x Nch'^^ to the data 
yielding the values k = 0.597 ± 0.009 ± 0.057 fm and k = 0.612 ± 0.007 ± 0.068 fm at 0.9 and 
at 7.0 TeV respectively. As can be seen, the two fitted curves almost coincide so that one can 




Figure 2: The i?iD values measured by the CMS collaboration [9] as a function of the charged-particle 
multiphcity Nch at y/s = 0.9 and 7.0 TeV . The dotted and sohd hnes represent the fit results of the 

1/3 

relation i?iD(-^ch) = kx Nj^ to the 0.9 and 7.0 TeV data respectively. The inner error bars represent 
the statistical uncertainties and the outer error bars the statistical and systematic uncertainties added 
in quadrature. 

safely conclude that -Rid is essentially independent of the energy but is correlated to the 
increase in the charge multiplicity as has already been observed in the studies of pp and light 
ion collisions [12] as well as in the e"'"e^ — )■ — hadrons annihilations [13]. 



3 The average charge multiphcity as a function of 
energy 

The average final state charge multiplicity (Ach) produced in pp and pp collisions is known to 
increase with energy as shown in the left part of Fig. 3 which was taken from Ref. [14]. The 
(A"ch) values shown were obtained from the charged tracks emitted in pp and pp collisions 
summed up over the entire phase space which was experimentally quite accessible in reactions 
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Figure 3: Left: A compilation, taken from [14], of the average charge multipUcity (A'ch) produced in 
pp and pp reactions versus ^/s covering the full phase space and within |?7| <1.5; Right: A compilation, 
taken from [15] of dNch/drj, versus y/s obtained from from pp and pp inelastic (INEL) reactions within 
1 7/1 < 0.5 and non— single diffractive collisions (NSD) within \ri\ < 1. The lines shown in the figures 
represent the fits to a power series in Ins. 

at a relatively low ^/s regions. However, it became more difficult as the energy increased with 
the hadron colliders to the values near 2 TeV at the pp Tevatron accelerator and even higher 
to the current value of 7 TeV supplied by the LHC at CERN where the analysis was restricted 
to a limited phase space defined by a range in the pseudorapidity rj. In this limited phase 
space one has still the possibility to determine dNch/drj i.e., the differential charge multiplicity 
density as a function of t]. To represent a dNch/dr] distribution one often uses its specific value 
at 7] = that is dNch/dri\r,=o- The right part of Fig. 3, taken from reference [15], shows the 
charge multiplicity density at ri=0 obtained in pp and pp collisions as a function of ^/s. A 
clear rise of both {Nch) and dNch/dri[ri=o with a/s is seen which can be parametrized by a 
power series in ln(s) as indicated in Fig. 3. Further properties concerning the charge 
multiplicity behavior as a function of energy can be found e.g. in reference [16]. Taking the 
(Nch) dependence on energy in GeV units [14] to be 3.105 x s^ i'^s and for Rid the one given 
by Eq. (2) one obtains 

i?iD ^ 0.24 + 0.392 X ln((A/'eh)) fm. (3) 

Combining the relation dN^Y,/ dr]\n=o = 2.716 - 0.307 x ln(s) + 0.0267 x In^(s) given in [15], 
where s is in GeV^ units, and the -Rid energy dependence given by Eq. (2) one has 

Rm ^ 1.08 + 0.43 x ^dNch/dri\^=Q - 1.83 fm. (4) 

which is valid from about dNch/ dr]\rj=o = 2. The Rid dependence on (Nch) and on 

dNch_/ dri\ri=Q are shown in Fig. 4. From this figure one may of course also estimate from the 

measured -Rid the values of {Nch) and dNch/ dri\r,=o at the same pp energy. 
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Figure 4: Left: The Rid dependence on the average charge multiphcity (A^ch) according to Eq. (3); 
Right: The Rid dependence on dNch/dri\r]=o according to Eq. (4). 

4 The Rid dependence on the number of jets 

Even though Rid is seen to rise with (A^ch), the BEC studies of the hadronic decays at the 
LEP e+e~ coUider have shown that the fundamental reason for the Rid rise is the increase of 
the number of the hadronic jets, A^jet- An example of this feature is illustrated in Fig. 5 where 
Rid and A results obtained from BEC analyzes of the Z° hadronic decay carried out by the 
OPAL collaboration at LEP [13] are plotted. In the left part of the figure Rid is seen to 
increase with A^ch while A decreases with it. In the right part of this figure one observes that 
the ii!iD values are consistently higher for 3-jet events than those derived from the 2-jet events 
while both sets of the -Rid values are only slightly dependent on A^ch- At the same time the A 
values do not reveal any obvious separation between the 3 and 2-jet events nor do they show a 
distinct dependence on the charge multiplicity. A similar result has also been obtained by the 
L3 Collaboration at LEP [17]. That the Rid value is indeed linked to the number of hadronic 
jets has already been successfully demonstrated in the study of the measured dependence of 
Rid on y/s in e'^e^ annihilations into hadrons [18]. 

5 The mean number of jets deduced from Rid and A 

A compilation of Rid values obtained from BEC analyzes of identical pions produced in heavy 
ion collisions [19] has shown that to a good approximation the relation Rid — 1-2 x A^^^ fm 
holds where A is the atomic number of the colliding nucleus. This it!iD dependence on A is 
consistent with the dependence of the radius of nuclei on A deduced from other types of 
experiments such as the electron nucleus scattering. The increase in Rid is attributed first to 
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Figure 5: The dependence of the chaoticity parameter A and the source dimension Rm on the observed 
charge multiphcity A?chand number of jets in the hadronic decays given in reference [13]. Left: A 
and Rid as a function of A^ch- Right: The dependence of A and Rid on the observed AT^jifor two jet 
(sohd points) and three jet events (open points). The straight lines in the figures represent Unear fits 
to the data points. 



the rise in the number of the hadron sources, that is the number of participating 
nucleon-nucleon coUisions, as A increases. The second feature which causes the Rid to rise is 
the absence of a BEC effect (i.e. A=0) in identical boson pair correlation emerging from 
different sources. This last feature has been experimentally verified e.g. in the process 
e'^e~ — > W'^W~ — > pions where the BEC effect is present in identical pion pairs emerging 
from the same W gauge boson but is absent when the pair was emitted from two different W 
gauge bosons [20] . Prom this follows that whenever identical bosons are collected for a BEC 
analysis from more than one source the over all correlation decreases which results, as briefiy 
elucidated below, in an increase in Rid ■ 

Inasmuch that the emitting sources are similar the Rid values obtained from BEC of two 
identical bosons emitted from and S^n sources are related by [3,21,22] 

Rm — 'Y'JT^^k , (5) 

where Dk and Dm are the reduction factors corresponding respectively to the Sk and Sm 
sources. It has further been shown that in the case that the sources have about the same 
charge multiplicity then for A'ch ^ 1 the reduction factor is simply given by 



This inverse dependence of Dj on the number of sources holds also when the occurrence of 
slight deficiencies in the multiplicity measurement are accounted for while the usage of the 
full dilution value implies that the system of sources fully overlap in the momentum 
space [3,23]. Thus the BEC deduced i?iD values from the Sk and sources are related by 
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Figure 6: The average number of hadronic jets produced in pp collisions multiplied by the BEC 
deduced A values as a function of ^/s in TeV units obtained via Eq. (9) relative to the average number 
of jets produced at ^ = 0.01 TeV multiplied by its A value. 

Sk Sk Xm Rk Sk Rkl Ajt 

As has been apphed successfully in the past [21] we also here, in particle reactions, associate 
the hadronic jets with the boson sources. As a consequence Eq. (7) can be rewritten as 



Thus from i?iD and A measurements via BEC analyzes at two different "p-p energies one should 
be able to estimate the ratio of the average number of hadronic jets (A'jet) produced over the 
full phase space at these two energies. Here it is worthwhile to note the the BEC method for 
jet analysis does not depend on a a particular choice of a jet finding algorithm and its 
variables. 

Next we explore the possibility to evaluate the energy dependence of the average number of 
hadronic jets produced in pp collisions. To this end, the energy dependence of i?iD given by 



Table 1: The average number of jets produced in pp collisions at the LHC energies relative to its 
value at 0.9 TeV estimated via Eqs (8) and (9). 



Ref. 


(TeV) 


Rib (fm) 
Measured 


A 

Measured 


(7Vj,,(^))/(7Vj,,(0.9)) 


Using Eq. (8) 


Using Eq. (9) 


[9] 


0.90 


1.56± 0.12 


0.62 ± 0.03 


1 


1 


[8] 


2.36 


1.99± 0.18 


0.66 ± 0.09 


1.20±0.22 


1.02±0.15 


[9] 


7.00 


1.89 ± 0.19 


0.62 ± 0.04 


1.21±0.18 


1.18±0.09 




14.0 


[2.01 ± 0.12]* 


[0.62 ± 0.03]** 




1.22±0.08 



*Extrapolated via Eq. (2) **Assumed 



Eq. (2) is incorporated into Eq. (8) with the result that 

Fig. 6 shows [Am x (A^jct(w^)) ]/[Afc x {Njct{k)) ] as a function of ^/s for pp collisions using in 
Eq. (9) with the values a = 1.64 fm and b = 0.14 fm and choosing (A'jet(/c)) to be the 
average number of hadronic jets expected at ^/s = 0.01 TeV. 

In Table 1 we summarize the Rid and A values recently measured in pp collisions at 0.9, 2.36 
and 7 TeV and their estimated values at 14 TeV, the originally planned LHC energy. Also 
given in the Table are the expected (N^ct) values relative to its value at 0.9 TeV. As can be 
seen, the results obtained for the average number of jets via Eqs (8) and (9) are consistent 
within errors apart from the value obtained at ^/s = 2.36 TeV where they differ considerable 
which is also the case for its Rib value which is seen in Fig. 1 to be on the higher side of the 
fitted one. Another feature to notice is the relatively slow increase of (N^et)- Namely, in an 
energy change by a factor of ~15, between 0.9 TeV to 14 TeV, the average number of 
hadronic jets increases just by about 22%. 



6 Summary and conclusions 

The one dimensional Rid deduced from Bose-Einstcin correlation of identical pions produced 
in pp collisions is seen to increase with the center of mass energy of the interacting baryons. 
This behavior can be parametrized by Rid = a + b x In(y^) fm where s is given in TeV^ units 
and with the fitted values of a = 1.64 fm and 6 = 0.14 fm. The Rid dependence on energy is 
then used to derive expressions for its dependence on the average charge multiplicity {Nch) 
and charge multiplicity density dNch/dri\,^=Q. It is further argued that the fundamental cause 
for the rise of Rid with ^/s is the increase with energy of the average number of outgoing 
hadronic jets produced in pp collisions. Two method are outlined for the estimation of the 
relative average true number of hadronic jets via BEC measurements of Rid and A which are 
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applied to the recently investigated pp collisions at the Large Hadron CoUider at CERN. An 
extrapolated value for the relative (A'"jet) value for the planned LHC pp collisions at 14 TeV is 
also given. 
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